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Abstract 

We present the Fortran code SDECAY, which calculates the decay widths and branch- 
ing ratios of all the supersymmetric particles in the Minimal Supersymmetric Stan- 
dard Model, including higher order effects. Besides the usual two-body decays of 
fermions and gauginos and the three-body decays of charginos, neutralinos and 
gluinos, we have also implemented the three-body decays of top squarks, and even 
the four-body decays of the top squark; the important loop-induced decay modes 
are also included. The QCD corrections to the two-body decays involving strongly 
interacting particles and the dominant components of the electroweak corrections 
to all decay modes are implemented. 



*The code can be obtained at the url: http://people.web.psi.ch/muehlleitner/SDECAY 
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1. Introduction 



The search for the new particles predicted by supersymmetric theories is a major goal of 
present and future colliders. A large theoretical effort has been devoted in the last two 
decades to determine the basic properties of these particles, as well as their decay modes 
and production mechanisms in collider experiments. Most of these studies have been 
performed in the framework of the Minimal Supersymmetric Standard Model (MSSM) 
[1-3], which has the minimal gauge structure, particle content and minimal interactions, 
and leads to a stable lightest supersymmetric particle (LSP). 

As is well known, it is a very complicated task to deal in an exhaustive way with all 
the properties of these new particles and to make detailed and complete phenomenological 
analyses and comparisons with the outcome of or expectations from experiments. This 
is mainly due to the fact that, even in the MSSM, there are more than a hundred new 
parameters in the most general case, and even if one constrains the model to have a 
viable phenomenology, there are still over 20 free parameters left to cope with. This large 
number of inputs enters in the evaluation of the masses of O(30) supersymmetric (SUSY) 
particles and Higgs bosons as well as their complicated couplings, which involve several 
non-trivial aspects, such as the mixing between different states, the Majorana nature of 
some particles and, if the aim is to be quite precise, the higher order corrections. One 
then has to calculate, in the most accurate way, the rates for the many possible decay 
modes and production processes at the various possible machines. 

Once SUSY particles are found, their properties are expected to be determined with 
an accuracy of a few per cent at the LHC [4] and a precision at the per cent level or 
below at future e + e~ linear colliders [5]. To match this expected experimental accuracy, 
we need to calculate the mass spectra, the various couplings, the decay branching ratios 
and the production cross sections with a rather high precision, i.e. including the higher 
order effects. Leaving aside the production processes which are dealt with by Monte Carlo 
event generators, and for which the cross sections have been calculated at next-to-leading 
order (NLO) in some cases [6], one therefore needs to achieve the following goals: 

- The physical (pole) masses and the various couplings [as well as the soft SUSY- 
breaking parameters which enter those] of the SUSY particles and the MSSM Higgs bosons 
need to be calculated very accurately. They must include all relevant features such as 
the masses of third-generation fermions, the mixing between the various states, and the 
radiative corrections when important. In constrained MSSMs two additional features 
need to be handled carefully: the renormalization group evolution (RGE) of parameters 
between the low-energy scale and the high-energy scale and the consistent implementation 
of radiative electroweak symmetry breaking (EWSB), i.e. the loop corrections to the 
effective scalar potential. There are several available RGE codes [7-10], which calculate 
the supersymmetric particle and Higgs boson masses as well as the soft SUSY-breaking 
parameters in the unconstrained and constrained MSSMs; they can be straightforwardly 
extended to allow for the calculation of the various and numerous couplings. 
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- All the possible two-body decay modes that can occur at the tree level [11] should 
be taken into account. These consist not only of decays of the inos 1 into fermion/sfermion 
pairs or sfermion decays into fermion/ino pairs, but also decays involving Higgs and gauge 
bosons in the final state. The QCD corrections, which are known to be rather large [12- 
14], need to be incorporated in the processes involving strongly interacting particles. In 
addition, some electroweak radiative corrections can be as large as the QCD corrections, 
and in principle, they should be taken into account. However, it is a tremendous task to 
calculate all these corrections, and they are available in the one-loop approximation [15] 
only for a very limited number of processes. Nevertheless, the bulk of these corrections, 
those stemming from the running of the gauge and Yukawa couplings and the running of 
some soft SUSY-breaking parameters, is available in the literature and can be incorporated 
with a minimum of effort. 

- All possibly important higher order decay modes must be included. They consist 
of the three-body decay modes of the charginos, neutralinos and the gluino into lighter 
inos and two massless fermions [16, 17], which are known to be important, but also of 
the three-body decays of third-generation sfermions, which have been shown recently to 
be possibly important in some kinematical and parameter configurations [18, 19]. This is 
particularly the case for the top squark, when its mass is smaller than the lighter chargino 
mass and the sum of the masses of the lightest neutralino and the top quark. In fact, 
even the four-body decay mode of the top squark [20] has been shown to possibly compete 
with the flavour changing neutral current (FCNC) and loop-induced decay into a charm 
quark and the lightest neutralino [21], which also has to be included. Another set of 
loop-induced decay modes which might play a prominent role are the radiative decays of 
the next-to-lightest neutralino into the lightest neutralino and a photon [22, 23] and, to 
a lesser extent, the decay of the gluino into a gluon and the LSP [24]. 

The Fortran code SDECAY, which we present in this report, deals with the decays of 
SUSY particles in the framework of the MSSM, and includes the most important higher 
order effects. It uses the RGE program SuSpect [8] for the calculation of the mass 
spectrum and the soft SUSY-breaking parameters [but of course, the program can be 
easily linked to any other RGE code] and evaluates the various couplings of the SUSY 
particles and MSSM Higgs bosons. It calculates the decay widths and the branching 
ratios of all the two-body decay modes, including the QCD corrections to the processes 
involving squarks and gluinos and the dominant electroweak effects to all processes. It also 
calculates the loop-induced two-body decay channels, as well as all the possibly important 
higher order decay modes: the three-body decays of charginos, neutralinos, gluinos and 
third generation sfermions and the four-body decays of the top squark 2 . 

^or simplicity, we will collectively call inos the charginos and neutralinos, and sometimes the gluinos. 

2 There are other programs which calculate the decay branching ratios of SUSY particles [7, 10, 25], 
and that also include some higher order effects. For instance, the loop decays and the three-body ino 
decays are included in ISASUGRA, SPHENO and SUSYGEN; the program SPHENO also deals with some three- 
body decays of the top squark. None of these programs includes the QCD corrections and the full set of 
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The program also calculates the decay widths and branching ratios of the heavy top 
quark. Besides the standard decay into a W boson and a bottom quark, the top quark 
decay widths into a charged Higgs boson and a bottom quark and into a top squark and a 
neutralino, are evaluated. The one-loop SUSY-QCD radiative corrections that are known 
[26] will be included in an upgraded version of the program. 

Besides the SuSpect files needed for the evaluation of the spectrum, the code contains 
only one source file sdecay.f, written in Fortran77, and one input file, sdecay.in, from 
which any choice of approximation in the calculation is driven [including or not the higher 
order corrections and/or decays, the choice of the various scales, the order of perturbation 
at which some couplings are calculated, etc.]. All results for the total decay widths and 
branching ratios are given in the output file sdecay . out, either in a simple form or in the 
SUSY Les Houches Accord (SLHA) [27] form [this choice can be made in the input file]. 
The program is very user-friendly, self-contained and it can easily be linked with other 
codes or Monte Carlo event generators. It is rather fast and flexible, thus allowing scans 
of the parameter space with several possible options and choices for model assumptions 
and approximations. 

The program SDECAY is of the same level of sophistication as the program HDECAY 
[28], which calculates the decay widths and the branching ratios of the [Standard Model 
and] MSSM Higgs bosons, including all relevant higher order effects. In fact, the three 
programs SuSpect, HDECAY and SDECAY, have many common features and subroutines, and 
are organized in a similar way. They provide a coherent, consistent and comprehensive 
description of the properties of the supersymmetric and Higgs particles in the MSSM, 
prior to the level of production which, as was mentioned previously, is the domain or 
"chasse gardee" of the Monte Carlo event generators. A light version, which combines 
these three programs, and which can be easily linked to any Monte Carlo generator, is 
under development and will appear quite soon [29]. 

The rest of this report is organized as follows. In the next section, we will summarize 
the main features of the MSSM that we will deal with, concerning the sparticle and Higgs 
boson masses and couplings, and the notation that we use in the program. In section 3, 
we discuss all the decay modes that are implemented and the way the higher order decays 
and the radiative corrections are included. In section 4, we discuss the main features and 
the structure of the program, briefly summarize how it works, and display the content of 
the input and output files. In section 5, a short conclusion will be given. 

2. The implementation of the MSSM 

The program SDECAY deals with the MSSM, i.e. with the basic assumptions of a: 

- minimal gauge group, the Standard Model SU(3)c x SU(2) L x U(l)y one, 
higher order decays of third-generation sfermions, though. 
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- minimal particle content: three generations of "chiral" sfermions f l LR [no right- 
handed sneutrinos] and two doublets of Higgs fields Hi and H 2 , 

- minimal set of couplings imposed by R-parity conservation to enforce baryon and 
lepton number conservation in a simple way and which leads to a stable LSP, 

- minimal set of soft SUSY-breaking parameters: gaugino mass terms M i: scalar mass 
terms and mf., a bilinear term B and trilinear sfermion couplings A{. 

For the superpotential and the minimal set of soft SUSY-breaking, which give all the 
interactions and couplings, we follow the notations that can be found in the users manual 
of the program SuSpect [8], which we use for the determination of the SUSY particle 
and Higgs boson spectra. To have a viable phenomenology and a reduced number of free 
parameters, we thus also make the following three assumptions: 

(i) All the soft SUSY-breaking parameters are real and therefore no new source of 
CP-violation is generated, in addition to the one from the CKM matrix. 

(ii) The matrices for the sfermion masses and for the trilinear couplings are all diagonal, 
implying the absence of FCNCs at the tree level. 

(Hi) The first and second sfermion generations are universal at low energy to cope with 
some severe experimental constraints [this is also motivated by the fact that we have 
neglected for simplicity all the masses of the first- and second-generation fermions 
which are small enough not to have any significant effect]. 

Making these three assumptions will lead to the so-called "phenomenological MSSM" 
(or pMSSM) discussed in [3], with 22 input parameters only: 

tan/5: the ratio of the vevs of the two-Higgs doublet fields. 

m 2 Hi ,m 2 H2 : the Higgs mass parameters squared. 

Mi, M 2 , M 3 : the bino, wino and gluino mass parameters. 

rriq, rriu R , m d - , mi, me R : the first /second-generation sfermion mass parameters. 

mQ,mi R ,mi R ,mi,mf R \ the third-generation sfermion mass parameters. 

A u ,A d ,A e : the first /second-generation trilinear couplings. 

A t , Ab, A T : the third-generation trilinear couplings. 

If one requires a proper electroweak symmetry breaking, the Higgs-higgsino (super- 
symmetric) mass parameter (up to a sign) and the soft SUSY-breaking bilinear Higgs 
term B are determined, given the above parameters [alternatively, one can trade the val- 
ues of m 2 Hl and m 2 H2 with the "more physical" pseudoscalar Higgs boson mass Ma and 
parameter \i, a possibility provided by the program SuSpect, which deals with all the 
aspects of EWSB]. 

In constrained models, such as minimal Supergravity (mSUGRA) [30], the gauge- 
mediated SUSY breaking (GMSB) [31] and anomaly-mediated SUSY breaking (AMSB) 
[32] models, most of the 22 soft SUSY-breaking input parameters of the pMSSM listed 
above are derived from a set of universal boundary conditions at the high-energy scale [the 
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Grand Unification (GUT) scale for the mSUGRA and AMSB models and the messenger 
scale for the GMSB model]. In the mSUGRA case, for instance, the entire set of soft 
SUSY-breaking parameters [and thus the superparticle and Higgs spectrum] is determined 
by the values of only five free parameters: a common gaugino mass TO1/2 = Mi, a universal 
scalar mass itlq — mj — m#. and a universal trilinear coupling Aq = Ai at the GUT scale, 
the sign of the higgsino parameter \x and tan f3. 

The low-energy soft SUSY-breaking parameters are then derived from the high-energy 
ones above through Renormalization Group Equations. [Note that the values of and 
Ma are obtained by requiring proper EWSB, which should be implemented, and the 
value of M G ut is defined as the scale where the three gauge coupling constants of the 
Standard Model (SM) unify]. One then proceeds to calculate the pole masses of the 
Higgs bosons and all the supersymmetric particles, including the possible mixing between 
the current states and the radiative corrections [up to two loops in some cases] when they 
are important. 

All these steps are performed by the program SuSpect and, for completeness, we repro- 
duce in Fig. 1 the general algorithm that is used in the code [8] . This iterative algorithm 
includes the various important steps of the calculation: the choice of SM input parameters 
at low energy [the gauge coupling constants and the pole masses of the third-generation 
fermions], the calculation of the running couplings including radiative corrections in the 
modified Dimensional Reduction DR scheme [which preserves SUSY] and their RG run- 
ning back and forth between the low and high scales, with the possibility of imposing the 
unification of the gauge couplings and the inclusion of SUSY thresholds in some cases, 
the RG evolution of the soft SUSY-breaking parameters from the high scale to the EWSB 
scale, the minimization of the one-loop effective potential and the determination of some 
important parameters, and finally the calculation of the particle masses including the 
diagonalization of the mass matrices and the radiative corrections. 

To be more specific, we provide below a summary list of the higher order effects that 
have been included in the program: 

- For the DR gauge and third-generation fermion Yukawa couplings, defined at the 
scale Mz, the full set of standard and SUSY corrections has been implemented according 
to the approach of Pierce, Bagger, Matchev and Zhang (PBMZ) [33]. There are two 
exceptions: in the case of sin 2 9w, the small SUSY particle contributions to the box dia- 
grams have been omitted, and in the case of m& and m T , only the QCD and the leading 
electroweak corrections at zero-momentum transfer have been incorporated [which, ac- 
cording to PBMZ, is a very good approximation]. In some cases, again according to the 
PBMZ approach, some important two-loop corrections have also been taken into account. 
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Choice of low energy input: a(Mz), sin 2 9w, &s(Mz), ^ttr i tan/3(M^) 
Radiative corrections g^ 3 (M z ), A^(M Z ), A^(M Z ), Ap(M z ) 
First iteration: no SUSY radiative corrections. 
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Two-loop RGE for a??, and A??, with choice: 91 92 ' V 3 / 5 

P 2/1,2,3 r,M M GUT ~ 2 x 10 16 GeV 

Include all SUSY thresholds via step functions in (3 functions. 

First iteration: unique threshold guessed. 



+ 

Choice of SUSY-breaking model (mSUGRA, GMSB, AMSB, or pMSSM). 
Fix your high-energy input (mSUGRA: m ,m 1 / 2 , A , sign(/x), etc.). 

t 



Run down with RGE to: f ° r f D 

-Mewsb for rrii, Mi, A iy //, B 

First iteration: guess for Mewsb = Mz- 



t 



H 2 , [iB — F non _i ineai (m Hl , rn H2 , tan (3, V\ oop ) 
Vi 00 p = Effective potential at one loop with all masses. 
First iteration: V\ oov not included 



Check of consistent EWSB (fi convergence, no tachyons, simple CCB/UFB, etc.) 

+ 



Diagonalization of mass matrices and calculation of masses / couplings 
Radiative corrections to the physical Higgs, sfermion, gaugino masses. 
First iteration: no radiative corrections. 



t 



Check of a reasonable spectrum: 
no tachyonic masses (from RGE, EWSB or mix), good LSP, etc.. 
not too much fine-tuning and sophisticated CCB/UFB conditions, 
agreement with experiment: Ap, (g — 2), b — > 57. 



Figure 1: Iterative algorithm for the calculation of the SUSY particle spectrum in SuSpect 
from the choice of input parameters (first step) to the check of the spectrum (last step). 
The EWSB "small" iteration on // the RG/RC "long" iteration are performed until a 
satisfactory convergence is reached. 
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- The RGEs have been used at the two-loop level for the gauge and Yukawa cou- 
plings, as well as for the three gaugino mass parameters M ± , M 2 , M 3 and the two Higgs 
mass parameters rriH 1 , ttih 2 • For the other soft SUSY-breaking parameters [essentially the 
sfermion mass parameters mj and the trilinear couplings Af], only the one-loop RGEs 
have been used. The GUT scale Mgut can be either defined as the value at which the 
two gauge couplings g± and g 2 unify, or can be set by hand at Mgut ~2x 10 16 GeV. 

- The EWSB has been implemented through the tadpole method. The full one- 
loop standard and SUSY contributions to the tadpoles have been taken into account. 
The dominant two-loop corrections, those stemming from QCD and the third-generation 
fermion Yukawa couplings, have also been included. The EWSB scale has been chosen to 
be the geometric mean of the two top squark masses, M E wsb = ( m t 1 m t 2 ) 1 ^ 2 - 

- The soft SUSY-breaking parameters and the parameter /i that we obtain are all 
DR parameters defined at the scale Mewsb- Using these parameters, the chargino and 
neutralino mass matrices are diagonalized with real matrices U /V and Z, respectively, to 
obtain the tree-level physical ino masses. The third-generation sfermion mass matrices 
are also diagonalized to obtain the tree-level sfermion masses, the mixing angles 0$, 0~ b and 
0f being DR parameters defined at the scale Mewsb- 

- The radiative corrections to the sfermion masses are again included according to 
PBMZ [33], i.e. only the QCD corrections for the superpartners of light quarks [including 
the bottom squark] plus the full QCD and electroweak corrections to the two top squarks; 
the small electroweak radiative corrections to the slepton masses [which according to 
PBMZ are at the level of one per cent] have been neglected. The full one-loop QCD 
radiative corrections to the gluino mass are incorporated, while in the chargino/neutralino 
case, the radiative corrections to the masses are simply included in the gaugino and 
higgsino limits, which is a very good approximation according to PBMZ. 

- The calculation of the masses of the MSSM Higgs bosons and the mixing angle a 
in the CP-even sector can be made using four routines that are available on the market: 
Subhpole [34], HMSUSY [35], FeynHiggsFast [36] and Hmasses [37]. In view of the ap- 
proximations used in SuSpect, it is more appropriate to use the routine Hmasses, which 
has the same level of approximation, since it includes the full one-loop radiative correc- 
tions, and the dominant two-loop corrections at order a s Xj and Aj, A/ being the Yukawa 
couplings of the third-generation fermions, f — t,b and r. 

Using the gauge couplings, the third-generation fermion masses and the soft SUSY- 
breaking parameters discussed above, we then proceed in the program SDECAY to calculate 
all the couplings of the SUSY particles and the MSSM Higgs bosons. In most cases, we use 
the Feynman rules given by Haber et al. in Ref. [2]. These couplings are thus defined in 
the DR scheme and evaluated at the scale Mewsb- Therefore, they already include some 
radiative corrections, and care should be taken when they are used in one-loop-corrected 
amplitudes to avoid double counting, as will be discussed later. 
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3. Decays of the SUSY particles 
3.1 Two-body decays at the tree level 

The main decay modes of sfermions will be into their partner fermions and neutrali- 
nos, as well as into their isospin partner fermions and charginos 3 

h - fxj- (i) 

In the case of squarks, when they are heavier than the gluino, they can also decay into 
gluino-quark final states 

qi->qg ■ (2) 

If the mass splitting between two sfermions of the same generation is large enough, as can 
be the case of the third-generation (t,b) and {y T ,f) isodoublets, the heavier sfermion can 
decay into the lighter one plus a gauge boson V = W, Z or a Higgs boson $ = h,H, A, if ± 

u - f;v, (3) 

fi - jj*. (4) 

The heavier neutralinos and charginos will decay into the lighter chargino and neu- 
tralino states and gauge or Higgs bosons 

X* - XjV (5) 
X* Xj® (6) 

and, if enough phase space is available, into fermion-sfermion pairs 

X* -> ffj ■ (7) 

For gluinos, when they are heavier than squarks, their only relevant decay channel will 
be into quark plus squark final states: 

g ->m ■ ( 8 ) 

In the case of the GMSB model , the lightest SUSY particle is the gravitino G. The 
next-to-lightest SUSY particle (NLSP) can be either the lightest neutralino or the lightest 
sfermion, in general the f\ [in which case, the decay Xi ~^ ^ is allowed kinematically] . If 
the NLSP is a slepton, its only allowed decay is into a lepton and a gravitino, £ — > £G, 
with a branching ratio of 1. If the NLSP is the lightest neutralino, there are several 
possible decays, Xi — ► lG, ZG and G$ with $ = h, H, A. 

3 Here and in the following, we collectively denote by \i the charginos and neutralinos, and we discard 
the distinction between the two isospin sfermion partners. 
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The program SDECAY calculates the partial widths and the branching ratios of all 
these decays, including all possible combinations. It first checks if the decay is phase- 
space-allowed and then uses the two-body simple formulae available in the literature; see 
for instance Refs. [17, 19]. The masses involved in the phase space are all pole masses 
[i.e., in the case of the SUSY and Higgs particles, the one- loop renormalized masses given 
by SuSpect]. The masses of the third-generation fermions [the only ones in the program 
that are assumed to be non-zero] are the on-shell masses when they enter the phase space 
but the running DR masses defined at the scale Mewsb when they enter the various 
couplings. This is also the case of all the soft SUSY-breaking parameters and the third- 
generation sfermion mixing angles that enter the couplings. Nevertheless, we have left as 
an option the possibility for the QCD coupling constant and the b, t Yukawa couplings, 
to be evaluated at the scale of the decaying superparticle or any other scale. However, in 
this case, only the standard QCD corrections are included in the running, following the 
approach of Ref . [38] . 

3.2 QCD corrections to the two-body decays 

The one-loop QCD corrections have been incorporated to the two-body decay processes 
involving (s) quarks or gluinos in the initial or final state; more specifically, they have been 
included in the following processes 

Qi -»• QXj (9) 
Qi - Qj® (10) 
qi -> qg and g -> q t q . (11) 

We have used the formulae given in Refs. [12], [13] and [14] for, respectively, the 
processes of eqs. (9), (10) and (11). The corrections for the process Xi ~^ QiQ that can 
be adapted from the ones of the reverse process, eq. (9), and the corrections to the decay 
Qi — > qjV will be implemented in the next version of the program. The Passarino-Veltman 
one-, two- and three-point functions for the loop amplitudes have been implemented 
using the formulae given in Ref. [39]. For the three-body phase-space integrals, when 
an additional gluon is emitted in the final state, we use the analytical formulae given in 
Refs. [12-14] and which involve at most Spence functions. 

A few remarks are worth making at this stage: 

(i) All the corrections have been incorporated in the DR scheme, while in the previous 
references some corrections were calculated in the MS scheme. In the latter case the 
differences between the two schemes have to be corrected by additional counterterms. 
The results for the physical observables of course will be the same in the two schemes up 
to the calculated order, as it should be. 

(ii) All the masses of the particles involved in the processes, in particular those of the 
strongly interacting particles, which in principle need to be renormalized, are pole masses, 
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i.e. are defined on the mass shell. The self-energies of these particles are defined in such 
a way that the residues at the poles are unity. 

(in) Because the top and bottom quark masses and the stop and sbottom mixing 
angles [and thus the trilinear couplings A t and A^] that are obtained using the program 
SuSpect, already include some one-loop contributions, some care has to be taken when 
dealing with the renormalization of these parameters and their one-loop counterterms to 
avoid double counting [in fact only the divergent pieces have to be included]. 

(iv) In the case of decays involving gluinos, which are strong interaction decays already 
at the tree level, a special treatment is needed for a s (fiR), where hr is the renormalization 
scale. The heavy particles, top quarks, squarks and gluinos should be removed from the 
fi R evolution of a s and decoupled for momenta smaller than their masses. 

Finally, as mentioned previously, the bulk of the electroweak radiative corrections 
which is due to the running of the gauge and third-generation fermion Yukawa couplings, is 
already taken into account since these parameters, when appearing in the amplitudes, have 
been evaluated at the EWSB scale. The remaining corrections [including photon emission 
in the initial or final state] are of the order of the electromagnetic coupling constant and 
should lead to corrections at the level of a few per cent only. These corrections can thus 
be safely neglected in a first stage. 

3.3 Loop-induced decays 

a) Radiative decays of the next-to-lightest neutralino (and gluino) 

If the mass splitting between the next-to-lightest neutralino and the lightest neutralino 
is very small, the two-body neutralino decays discussed previously are not allowed kine- 
matically; the virtuality of the exchanged particles in the possible three-body decay modes 
[to be discussed later] is so large that the loop-induced decays of at least x° m to the LSP 
Xi and a photon [22, 23] 

xl - X\l (12) 

might be relevant. This decay is induced by triangle diagrams involving the contribution 
of virtual charginos together with W and charged Higgs bosons, and contributions with 
charged sfermion/fermion loops. It is of 0(a 3 ) in the electroweak coupling, but can have 
a sizeable branching ratio [with respect to the three-body decays] in some areas of the 
MSSM parameter space and in some corners of the phase space. 

For the full analytical formulae of the decay amplitudes, we use the ones that were 
given in Ref. [22], which are closer to our Lagrangian convention and leads to an easier 
implementation in our program. All diagrams and contributions have been taken into 
account. 

For completeness, and although this mode is never very important in the MSSM, we 
have extended the previous calculation to the case where the gluino decays into a gluon 
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and the lightest neutralino [24] 
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(13) 



which is mediated by only a subset of the diagrams involved in the previous decay mode 
[that is, only the squark-quark loop contributions]. 

b) Loop-induced decay of the lightest top squark 

The heaviness of the top quark leads to distinct phenomenological features for the 
decays of its scalar partners. Indeed, while the other squarks can decay directly into 
(almost) massless quarks and the lightest neutralino Xii which is always kinematically 
accessible, since in general, in the MSSM the neutralino Xi is assumed to be the LSP, 
the decay channels t« — > t\\ are kinematically closed for m^. < m t + m x o- If, in addition, 
m t x < iTT-b + m x +, the decay mode into a chargino and a b quark, ti — > bxf , is not accessible 



and the only two-body decay channel that would be allowed is the loop-induced and FCNC 
decay [21]: 



This mode is mediated by one-loop diagrams: vertex diagrams as well as squark and quark 
self-energy diagrams [where bottom squarks, charginos, charged W and Higgs bosons are 
running in the loops]. The flavour transition b — > c occurs through the charged currents. 
Adding the various contributions, a divergence is left out, which must be subtracted by 
adding a counterterm to the scalar self-mass diagrams. When working in mSUGRA-type 
models, where the squark masses are unified at the GUT scale, the divergence is subtracted 
using a soft counterterm at Aqut, generating a large residual logarithm log(AQ UT /M^) ~ 
65 in the amplitude. This logarithm gives the leading contribution to the ti — > cx? 
amplitude and makes the decay width rather large, although it is suppressed by the CKM 
matrix element V c b ~ 0.05 and the (running) b quark mass squared ml ~ (3 GeV) 2 . It 
is this approach that we have implemented in the program SDECAY and we have used the 
approximate formulae of Ref. [21] [the exact expressions of the loop amplitudes are not 
yet available]. 

However, there are scenarios in which the decay rate T(ti — > cxi) can be rather small: 
(i) First, the large logarithm log (Aq UT /M^) ~ 65 appears only because the choice of the 
renormalization condition is made at Aqut, but in a general MSSM, where the squark 
masses are not unified at some very high scale, one could chose a low-energy counterterm; 
in this case no large logarithm would appear, (ii) If the lightest top squark is a pure 
right-handed state [as favoured by the constraints from high-precision electroweak data], 
the amplitude involves only one component, which can be made small by choosing small 
values of the trilinear coupling Af, and/or large values of the (common) SUSY-breaking 
scalar mass m q . (Hi) Even in the presence of stop mixing, for a given choice of MSSM 
parameters, large cancellations can occur between the various terms in the loop amplitude; 
in addition, the t-ti~Xi coupling, which enters as a global factor, can be very tiny. 



ti 




(14) 
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Thus, the decay rate T(tx — > c%?) might be very small, opening the possibility for the 
three-body and even four-body decay modes, which will be discussed now, to dominate. 

3.4 Multibody decay modes 

If the tree-level two-body decay modes discussed previously are kinematically closed, i.e. 
equal to zero, multibody final state channels [as well as the loop decays which have been 
discussed previously] will be the dominant decays and should be considered. There are 
three varieties of such decays with which the program SDECAY deals, and which will be 
discussed below. In all cases, we make a clear separation between the two-body and 
multibody decays: that is, in the higher order decay modes, we do not include the total 
decay widths of the virtually exchanged (s)particles in their respective propagators, so as 
to make smooth transitions between the two-body and the multibody decay modes 4 . 

a) Three-body decays of charginos, neutralinos and gluinos 

When the two-body channels of the charginos and neutralinos given in eqs. (5-7) are 
kinematically closed, the particles can decay into three-body final states involving a lighter 
ino and two massless fermions: 

X^ - Xjff ■ (15) 

In the past, these decays have been discussed in particular for the lightest chargino and for 
the next-to-lightest neutralino, when decaying into the LSP neutralino and two fermions, 
Xi — > Xiff an d X2 ~* Xiff- I 11 our case, we will not assume that the initial inos are 
xf and Xi an d the final neutralino is the LSP Xii but any of the charginos xt = 1? 2) 
and neutralinos Xi — !>•••, 4) to cover all situations and the possibility of cascade 
decays. These decays proceed through gauge boson exchange [V = W and Z for xt an d 
Xi decays, respectively], Higgs boson exchange [$ = H + for xt decays and $ = H,h, A 
for x°i decays] and sfermion exchange in the t- and w-channels [the flavour is fixed by the 
sfermion-fermion and final neutralino vertex]. 

The gluino can also undergo three-body decays when the two-body decays of eq. (8) 
are kinematically forbidden: 

9 -»• XQQ , (16) 

and only the channels with t- and w-channel squark exchange will be present in this case; 
the partial widths can be straightforwardly derived from those of the neutralino decays, 

4 This aspect leads to very complicated technical problems and numerical instabilities. Indeed, once 
the total decay widths of the virtual particles are included, the phase-space integrals cannot be made 
analytically and one has to resort to a numerical integration which gives less precise results and heavily 
slows down the program. In addition, there are issues of gauge invariance that are not yet settled in this 
case, and which we did not want to address at the moment. 
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with the appropriate change of the couplings and the QCD factors. We also implemented 

the decay 

g - hbW (17) 

which can be important in regions of the parameter space where it is kinematically allowed 
[40]. Furthermore the decay [41] 

g -> hbH- (18) 

has been taken into account. 

In fact, the heavier charginos and neutralinos can also decay [in particular, in models 
with non-universal gaugino masses at high scale] into gluinos and quark-antiquark pairs 

Xi -> gqq , (19) 

and the amplitude can be adapted from the one of the reverse process eq. (16) discussed 
above. 

The program SDECAY calculates the partial widths and branching ratios of all these 
three-body decay channels, taking into account all the possible contributions of the vir- 
tual particles. The radiatively corrected Yukawa couplings of third-generation fermions, 
the mixing pattern for their sfermion partners and the masses of the sparticles and 
gauge/Higgs bosons involved in the processes are taken into account. In fact, even the 
masses of the fermion final states have been taken into account, since finite fermion masses 
are needed in some cases 5 . We have used the analytical formulae of Ref. [17] for the matrix 
elements squared and integrated numerically over the three-body phase space. 

b) Three-body decays of third-generation sfermions 

For relatively heavier top squarks, when their masses are larger than the mass of 
the Xi neutralino and the mass of the W boson, the boson and/or the sfermion /*, 
there will be the possibility of the three-body decay modes of the t±. These modes have 
been discussed in Refs. [18, 19]. When kinematically possible, some of these channels 
can dominate over the loop-induced U — > cx\ mode in rather large areas of the MSSM 
parameter space. 

For TOfj > M W (h±) + m x\i ^ ne three-body decay channels into W or charged Higgs 
bosons, a bottom quark and the LSP neutralino, 

U-*_ bW+xl , bH + x°i (20) 

5 The approximation of zero mass is rather good for all light fermion final states, except for 6-quarks 
and r-lcptons when x% an( ^ xt have masses close to the Xi mass; the zero-mass approximation would 
be very bad for top quark final states. Nevertheless, in mSUGRA-type models, if the three-body decays 
X2 ~ * Xitt an d Xi ~^ Xitb are kinematically allowed, they will not play a major role, since the charginos 
and neutralinos will have enough phase space to decay first into the two-body channels x\ ~ * Xi^, xV 1 
[and possibly x\H and x?A] and xt ~^ x\W [and possibly \\H + ], which will be largely dominating. 
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can be accessible and have been shown to be [at least for the one with W boson final 
states] often dominant in the case where < m t + m x o and m h + m x + . In addition, 
if sleptons are lighter than squarks [as is often the case in models with a common scalar 
mass at the GUT scale, such as the mSUGRA model] the modes 

U -> bl + ui and/or bl + vi (21) 

become possible. In the case of the lightest top squarks, they can be largely dominating 
over the loop-induced cxi mode. 

In fact, these three-body decay modes are important not only for the lightest top 
squark, but also for the heavier one. In addition, there is another possibility which is the 
decay of the top squarks into a fermion-antifermion pair and the lightest b state [which can 
become the lightest scalar quark in the case where tan/3 is very large], which is mediated 
by the virtual exchange of W and H + bosons: 

U - hff . (22) 

For the heavier top squark t 2 , another possibility would be the three-body decay into the 
lightest top squark and a fermion pair [with / ^ b] through the exchange of the Z and 
the MSSM neutral Higgs bosons [the CP-even h, H and the CP-odd A bosons], 

t 2 - hff- (23) 

These modes apply also for the charged decays of heavier bottom squarks into top squarks 
(and vice versa) which, as previously, occur through W and H + boson exchanges 

h - % ff ■ (24) 

[If the mass splitting between the initial and final scalar eigenstates is large enough, the 
gauge and Higgs bosons become real, and we have the two-body decays into gauge and 
Higgs bosons which have been discussed previously] 

For bb final states, one needs to include, in the case of t 2 — > t\ bb, the contributions 
of the exchange of the two charginos states Xi2- This is also the case of the decay mode 
b 2 — > bi bb, where one has, in addition, the virtual exchange of neutralinos and gluinos, 
which have to be taken into account. The latter process is a generalization [since the 
mixing pattern is more complicated] of the decay modes of first- and second-generation 
squarks into light scalar bottoms [19], and would be in competition with at least the 
two-body mode b 2 — > bxi- The latter channel is always open since Xi is the LSP, but the 
b-b 2 -Xi coupling can be small, leaving the possibility to the three-body mode to occur at 
a sizeable rate. 

SDECAY evaluates all the three-body decay modes of the top squarks discussed above, 
when the corresponding two-body decay channels are kinematically closed. In all cases, 
the analytical expressions for the Dalitz plot densities in terms of the energies of the 
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final fermions [19], which take account of all contributing channels, mixing and masses 
[including non-zero third-generation fermion masses], are integrated numerically. As in 
the case of the inos, the total decay widths of the exchanged particles are not included 
in the propagators of the virtual particles. The three-body decays of the bottom squarks 
will be included in the upgraded version of SDECAY. 

c) The four-body decay of the top squark 

If the previous three-body decay modes of the top squarks are kinematically not ac- 
cessible, the main t\ decay channel is then expected to be the loop-induced and flavor- 
changing decay into a charm quark and the LSP, t\ — > cx\- However, there is another 
decay mode that is possible in the MSSM, even if t\ is the next-to-lightest SUSY particle 
[provided that > m x o + m b ]: the four-body decay into a 6-quark, the LSP and two 
massless fermions [20]: 

h - b X °Jf' • (25) 

This mode is of the same order of perturbation theory as the decay t\ — > cxi, i-e. 0(a 3 ); 
in principle, it can therefore compete with the latter channel. 

The four-body decay mode proceeds through several diagrams: there are first the W- 
boson exchange diagrams with virtual t, b and xt,2 states, a similar set of diagrams is 
obtained by replacing the lU-boson by the if + -boson and a third type of diagrams consist 
of up and down type slepton and first /second-generation squark exchanges. The decay 
rate has been calculated in Ref. [20] taking into account all diagrams and interferences. 
The various contributions can be summarized as: 

i) Because in the MSSM, has a mass larger than M w and has tiny Yukawa cou- 
plings to light fermions, it does not give rise to large contributions. The squark exchange 
diagrams give also small contributions since squarks are expected to be much heavier 
than the t\ state. The contribution of the diagram with an exchanged t-quark is only 
important if the stop mass is of the order of m t + m x o > (9(250 GeV). 

ii) In contrast to squarks, slepton [and especially z>] exchange diagrams might give 
substantial contributions, since / masses of (9(100 GeV) are still experimentally allowed 
[42]. In fact, when the difference between and m x + and m^- is not large, this diagram 
will give the dominant contribution to the four-body decay mode. 

Hi) The most significant contributions to the four-body decay mode will come in gen- 
eral from the diagram in which the lightest chargino xt an d the W boson are exchanged, 
when the virtuality of the chargino is not too large. In particular, for an exchanged xt 
with a mass not much larger than 100 GeV, the decay width can be substantial even for 
top squark masses of the order of 100 GeV. 

Thus, a good approximation [especially for a light top squark mf ~ (9(100 GeV)] is 
to take into account only the top quark, the lightest chargino and the slepton exchange 
diagrams. All the other contributions, in particular the if ± contribution, can be safely 
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neglected in general. This is the approach that we choose in SDECAY: we include all the 
contributions and the interference effects, except for the small contribution which is 
very lengthy and time consuming to evaluate. The four-body phase-space integrals are 
evaluated using the program Rambo. 

3.5 Decays of the top quark 

The decays of the heavy top quark in the MSSM that have also been included in the 
program consist of the standard decay into a W boson and a bottom quark, 

t -> bW + (26) 

but also, if allowed by phase space, of the decays into a bottom quark and an H ± boson 
and into a lighter top squark and a neutralino 

t -> bH + and hxl (27) 

In addition, the one-loop QCD corrections [26] will be included in the next version of the 
program. 

4. Running SDECAY 

4.1 Basic facts about SDECAY 

Besides the files of the program SuSpect, i.e. the input file suspect2.in [where one can 
select the model to be investigated, the accuracy of the spectrum algorithm, the input 
data (SM fermion masses and gauge couplings) as well as SUSY and soft SUSY-breaking 
parameters] and the main Fortran routine suspect 2 . f [where some internally documented 
changes have been performed for a fully consistent calculation, and the calling routine 
suspect2_call .f is not needed], where the calculation of the spectrum is performed and 
which needs the routines subhJidec . f , feynhiggs.f, hmsusy.f and Hmasses for the 
calculation of the Higgs boson masses, the program SDECAY is composed of three files: 

1) The input file sdecay . in: in this file, one can choose the accuracy of the algorithm 
and make the choice of the various options: whether QCD corrections and multibody or 
loop decays are included, which scales are used for the couplings, the number of loops in 
their running and if top and GMSB decays are to be evaluated or not. 

2) The main routine sdecay. f where the couplings of the SUSY particles and Higgs 
bosons are evaluated and the decay branching ratios and total widths are calculated. This 
routine is self-contained and includes all the necessary files for the calculation. 

3) The output file sdecay. out: this file is generated at each run of the program and 
gives the results for the output branching ratios and total decay widths. Two formats are 
possible for this file: either the Higgs and SUSY particles are denoted in a simple and 
transparent form, or the PDG notation is used. In addition, the masses of the SUSY and 
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Higgs particles, the mixing matrices and the gauge and third-generation couplings at the 
EWSB or a chosen scale are given. [An additional output suspect 2 . out is also generated 
for the spectrum and also includes the soft SUSY-breaking terms]. 

The routine sdecay.f consists of about 30.000 lines of code and takes about 1 Mo of 
memory, while the input file only has a few dozen lines (most of them comments). The 
accompanying routines for the calculation of the sparticle and Higgs masses, which are 
provided separately, have in total a comparable size. The complete executable file takes 
about 2.5 Mo of disk space. The running time for a typical model point, for instance the 
mSUGRA point discussed below, is a few seconds on a PC with a 1 GHz processor. 

The Fortran files have to be compiled altogether and, running for instance on a PC, 
the compilation and link commands are [they are provided in a makefile]: 

OBJS = suspect2.o subh_hdec.o feynhiggs.o hmsusy.o sdecay.o 

FC=f77 

.f .o: 

$(FC) -c $*.f 
sdecay: $(0BJS) 

$(FC) $(0BJS) -o run 

Thus, the program SDECAY has the following structure: 

- It reads all the inputs in the files suspect2.in and sdecay. in. 

- It calculates the sparticle and Higgs boson masses as well as all the soft SUSY- 
breaking parameters using the program SuSpect. 

- It calls the subroutine common_ini where all parameters necessary for the calculation 
of the couplings and decay widths are set. 

- It calls the subroutine couplings where all couplings necessary for the calculation 
of the decay widths are evaluated. 

- It calls the subroutines for the two-body decays (2), the three-body decays (3) and 
loop decays (1) and for the stop four-body decay (4) calculation of the total widths 
and branching ratios of the respective decaying particle: 

chargino decays: subroutines: char2bod (2), xintegchipm (3) 

neutralino decays: subroutines: neut2bod (2), xintegneut (3), 



gluino decays: 



sup decays: 
sdown decays : 
stop decays: 



neutraddecay (1) 
subroutines : glui2bod (2) , xinteggo (3) , 

gluiraddecay (1) 
subroutine : sup2bod (2) 
subroutine : sdown2bod (2) 
subroutines: st2bod (2), xintegstop (3), 

hikasakobl (1) , st4bod (4) 
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sbottom decays: subroutine : sb2bod (2) 

selectron decays: subroutine : sel2bod (2) 

sneutrino_el decays: subroutine : snel2bod (2) 
stau decays: subroutine : stau2bod (2) 

sneutrino_tau decays: subroutine : sntau2bod (2) 

The routines call several help functions and subroutines for the loop decays, the 
QCD corrections, and some matrix elements for the multibody decays. They can 
be found at the end of the program with some comments specifying their purposes 
and their main features. 

- It writes in the output file sdecay.out where you can choose two versions: the 
output a la Les Houches Accord or one that is easier to read. 

In the next subsections, we will exhibit the input and output files, taking the example 
of an mSUGRA benchmark point from the Snowmass Points and Slopes [43], the so-called 
SPSla point, with the inputs at the high scale: 

m = 100 GeV , m 1/2 = 250 GeV , A = -100 , tan/5 = 10 , \i > . 

The input and output files are self-explanatory and will not be commented further. 

4.2 The input file 

SDECAY INPUT FILE 



* Choice of the output, Les Houches Accord (1) or simple (0): 
1 

* Include (1) or not (0) the QCD corrections to the 2-body decay widths: 


* Include (1) or not (0) the multi-body decays for inos and stops: 
1 

* Include (1) or not (0) the loop induced decays for the gluino, 
the neutralinos and stopl : 

1 

* Include (1) or not (0) the SUSY decays of the top quark: 
1 



* Include (1) or not (0) the possible decays of the NLSP in GMSB models: 
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(ichoice(l) has to be set 11 in suspect2 . in. ) 





* Scheme in which the running alphas and quark masses are calculated 
if the scale is not the scale of EWSB: 

(If QCD corrections are included, the DR_bar scheme has to be used.) 
(MS_bar scheme) and 1 (DR_bar scheme) . 



# Scale at which the scale dependent couplings are calculated: 

1: EWSB scale, 2: mass of the decaying sparticle, 3: user choice 

1 

# Scale of the scale dependent couplings if chosen by the user 
(in GeV) : 

100. DO 

# Number of loops for the calculation of the running couplings 
2 

4.3 Output file according to SLHA 

# SUSY Les Houches Accord - MSSM Spectrum + Decays 

# SDECAY 1.0 

# Authors: M.Muhlleitner, A.Djouadi and Y.Mambrini 

# In case of problems please send an email to 

# margarete.muehlleitner@psi.ch 

# djouadi@lpm.univ-montp2.fr 

# mambrini@delta.ft.uam.es 
# 

# If not stated otherwise all couplings and masses 

# are given at the scale of the electroweak symmetry 

# breaking Q= . 46296529E+03 
# 

# 

BLOCK MASS # Mass Spectrum 

# PDG code mass particle 



1 



25 



1 . 14365068E+02 



# h 



35 



3.91956602E+02 



# H 



36 



3.92191912E+02 



# A 



37 



4.00353329E+02 



# H+ 
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1000001 


5, 


.69828109E+02 


# 


~d_L 


2000001 


5, 


.43157826E+02 


# 


~d_R 


1000002 


5, 


.64244153E+02 


# 


~u_L 


2000002 


5, 


.44093303E+02 


# 


~u_R 


1000003 


5. 


.69828109E+02 


# 


~s_L 


2000003 


5, 


.43157826E+02 


# 


~s_R 


1000004 


5, 


.64244153E+02 


# 


~c_L 


2000004 


5, 


.44093303E+02 


# 


~c_R 


1000005 


5, 


. 16713072E+02 


# 


~b_l 


2000005 


5. 


.44166483E+02 


# 


~b_2 


1000006 


4. 


.00256829E+02 


# 


~t_l 


2000006 


5 


.80537860E+02 


# 


~t_2 


1000011 


2. 


.03724637E+02 


# 


~e_L 


2000011 


1 , 


.45386789E+02 


# 


~e_R 


1000012 


1. 


.87810224E+02 


# 


~nu_eL 


1000013 


2 


.03724637E+02 


# 


~mu_L 


2000013 


1. 


.45386789E+02 


# 


~mu_R 


1000014 


1 , 


.87810224E+02 


# 


~nu_muL 


1000015 


1. 


.36395255E+02 


# 


~tau_l 


2000015 


2. 


.07528018E+02 


# 


~tau_2 


1000016 


1. 


.86942145E+02 


# 


~nu_tauL 


1000021 


6 


.03561040E+02 


# 


~g 


1000022 


9. 


.89200644E+01 


# 


~chi_10 


1000023 


1. 


.76248916E+02 


# 


~chi_20 


1000025 


-3, 


.57870532E+02 


# 


~chi_30 


1000035 


3, 


.77017717E+02 


# 


~chi_40 


1000024 


1, 


.75568747E+02 


# 


~chi_l+ 


1000037 


3, 


.77194407E+02 


# 


~chi_2+ 



# 

BLOCK MIX # Neutralino Mixing Matrix 



1 


1 


9. 


. 84337446E- 


-01 


# 


N_ 


.11 


1 


2 


-6. 


25292026E- 


-02 


# 


N_ 


.12 


1 


3 


1. 


.54406991E- 


-01 


# 


N_ 


.13 


1 


4 


-5. 


.76920450E- 


-02 


# 


M_ 


.14 


2 


1 


1. 


. 12653886E- 


-01 


# 


N_ 


.21 


2 


2 


9. 


40622674E- 


-01 


# 


M_ 


.22 


2 


3 


-2. 


.77623463E- 


-01 


# 


N_ 


.23 


2 


4 


1, 


.59572240E- 


-01 


# 


M_ 


.24 


3 


1 


6. 


. 15283938E- 


-02 


# 


N_ 


.31 


3 


2 


-9. 


. 19478832E- 


-02 


# 


N_ 


.32 


3 


3 


-6. 


.94945017E- 


-01 


# 


N_ 


.33 
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3 


4 


-7. 


10500716E- 


-01 


# 


N_ 


.34 


4 


1 


1. 


. 20843498E- 


-01 


# 


M_ 


.41 


4 


2 


-3. 


.20725226E- 


-01 


# 


M_ 


.42 


4 


3 


-6. 


45085357E- 


-01 


# 


N_ 


.43 


4 


4 


6. 


.82932691E- 


-01 


# 


M_ 


.44 



# 

BLOCK UMIX # Chargino Mixing Matrix U 



1 


1 


-9. 


. 12748331E- 


-01 


# 


U_ll 


1 


2 


4. 


08522319E- 


-01 


# 


U_12 


2 


1 


4. 


.08522319E- 


-01 


# 


U_21 


2 


2 


9. 


. 12748331E- 


-01 


# 


U_22 


.OCK VMIX 


# Chargino 


Mixing 


Matrix V 


1 


1 


-9. 


71766502E- 


-01 


# 


V_ll 


1 


2 


2. 


. 35944624E- 


-01 


# 


V_12 


2 


1 


2. 


. 35944624E- 


-01 


# 


V_21 


2 


2 


9. 


71766502E- 


-01 


# 


V_22 



# 

BLOCK STOPMIX # Stop Mixing Matrix 

1 1 5.50903293E-01 # cos(theta_t) 

1 2 8.34569087E-01 # sin(theta.t) 

2 1 -8.34569087E-01 # -sin(theta.t) 
2 2 5.50903293E-01 # cos(theta_t) 

# 

BLOCK SBOTMIX # Sbottom Mixing Matrix 
1 1 9.21378487E-01 # cos(theta.b) 

1 2 3.88666546E-01 # sin(theta.b) 

2 1 -3.88666546E-01 # -sin(theta.b) 
2 2 9.21378487E-01 # cos(theta.b) 

# 

BLOCK STAUMIX # Stau Mixing Matrix 

1 1 2.80106527E-01 # cos (theta.tau) 

1 2 9.59968923E-01 # sin(theta.tau) 

2 1 -9.59968923E-01 # -sin(theta.tau) 
2 2 2.80106527E-01 # cos (theta.tau) 

# 

BLOCK ALPHA # Higgs mixing 

-1 . 13249720E-01 # Mixing angle in the neutral Higgs boson sector 

# 

BLOCK HMIX 0= 4.62965294E+02 # DRbar Higgs Mixing Parameters 
1 3.51486069E+02 # mu 
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BLOCK GAUGE Q= 4 . 62965294E+02 # The gauge couplings 

1 3.62163400E-01 # gprime(Q) DRbar 

2 6.46905504E-01 # g(Q) DRbar 

3 1.09847635E+00 # g3(Q) DRbar 

# 

BLOCK AU, AD, AE Q= 4 . 62965294E+02 # The trilinear couplings 
3 3 -5.11225438E+02 # A_t DRbar 
3 3 -7.92584896E+02 # A_b DRbar 
3 3 -2.54143182E+02 # A.tau DRbar 

# 

BLOCK Y_X,A_X Q= 4 . 62965294E+02 # The Yukawa couplings 
3 3 1.35709601E+00 # y_t DRbar 
3 3 2.10275093E-01 # y_b DRbar 
3 3 1.55612239E-01 # y_tau DRbar 

# 
# 
# 

# ================= 

# | The decay table I 

# ================= 

# 

# The multi-body decays for the inos and sfermions are included. 
# 

# The loop induced decays for the gluino, neutralinos and stops 

# are included. 
# 

# The SUSY decays of the top quark are included. 
# 

# 

# PDG Width 

DECAY 6 1.50609870E+00 # top decays 





BR 


NDA 


ID1 


ID2 














1. 


.00000000E+00 


2 


5 


24 


# 


BR(t 


-> 


b 




W+) 


0, 


•00000000E+00 


2 


5 


37 


# 


BR(t 


-> 


b 




H+) 


0. 


•00000000E+00 


2 


1000006 


1000022 


# 


BR(t 


-> 


~t_ 


.1 


~chi_10) 


0, 


.00000000E+00 


2 


1000006 


1000023 


# 


BR(t 


-> 


~t_ 


.1 


~chi_20) 


0, 


.00000000E+00 


2 


1000006 


1000025 


# 


BR(t 


-> 


~t_ 


.1 


~chi_30) 


0, 


.00000000E+00 


2 


1000006 


1000035 


# 


BR(t 


-> 


~t_ 


.1 


~chi_40) 


0, 


.00000000E+00 


2 


2000006 


1000022 


# 


BR(t 


-> 


~t_ 


_2 


~chi_10) 


0, 


.00000000E+00 


2 


2000006 


1000023 


# 


BR(t 


-> 


~t. 


_2 


~chi_20) 
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O.OOOOOOOOE+00 
. OOOOOOOOE+00 

# 

# PDG 
DECAY 1000021 

# BR 
1.76180098E-02 
1.76180098E-02 
5.39508838E-02 
5.39508838E-02 
2.37062902E-02 
2.37062902E-02 
5.23780813E-02 
5.23780813E-02 
1.76180098E-02 
1.76180098E-02 
5.39508838E-02 
5.39508838E-02 
2.37062902E-02 
2.37062902E-02 
5.23780813E-02 
5.23780813E-02 
1.01674937E-01 
1.01674937E-01 
5.53319969E-02 
5.53319969E-02 
4.76865361E-02 
4.76865361E-02 
. 00000000E+00 
. 00000000E+00 

# 

# PDG 
DECAY 1000006 

# BR 
1.76720929E-01 
1.26493082E-01 
. 00000000E+00 
. 00000000E+00 
6.78820856E-01 
1.79651325E-02 
0.00000000E+00 



2 2000006 
2 2000006 



Width 
4.85459975E+00 



NDA 


ID1 


2 


1000001 


2 


-1000001 


2 


2000001 


2 


-2000001 


2 


1000002 


2 


-1000002 


2 


2000002 


2 


-2000002 


2 


1000003 


2 


-1000003 


2 


2000003 


2 


-2000003 


2 


1000004 


2 


-1000004 


2 


2000004 


2 


-2000004 


2 


1000005 


2 


-1000005 


2 


2000005 


2 


-2000005 


2 


1000006 


2 


-1000006 


2 


2000006 


2 


-2000006 



Width 



2.06530974E+00 



NDA 


ID1 


2 


1000022 


2 


1000023 


2 


1000025 


2 


1000035 


2 


1000024 


2 


1000037 


2 


1000021 



1000025 # BR(t -> ~t_2 ~chi_30) 
1000035 # BR(t -> ~t_2 ~chi_40) 



# gluino decays 
ID2 



-1 


# 


BR( 


~g 


-> 


~d_ 


.L 


db) 


1 


# 


BR( 


~g 


-> 


~d_ 


.L* 


d ) 


-1 


# 


BR( 


~g 


-> 


~d_ 


.R 


db) 


1 


# 


BR( 


~g 


-> 


~d_ 


_R* 


d ) 


-2 


# 


BR( 


~g 


-> 


~u_ 


_L 


ub) 


2 


# 


BR( 


~g 


-> 


~u_ 


_L* 


u ) 


-2 


# 


BR( 


~g 


-> 


~u. 


_R 


ub) 


2 


# 


BR( 


~g 


-> 


~u. 


_R* 


u ) 


-3 


# 


BR( 


~g 


-> 


~s_ 


.L 


sb) 


3 


# 


BR( 


~g 


-> 


~s. 


_L* 


s ) 


-3 


# 


BR( 


~g 


-> 


~s_ 


_R 


sb) 


3 


# 


BR( 


~g 


-> 


~s. 


_R* 


s ) 


-4 


# 


BR( 


~g 


-> 


~c_ 


.L 


cb) 


4 


# 


BR( 


~g 


-> 


~c_ 


.L* 


c ) 


-4 


# 


BR( 


~g 


-> 


~c_ 


.R 


cb) 


4 


# 


BR( 


~g 


-> 


~c_ 


_R* 


c ) 


-5 


# 


BR( 


~g 


-> 


~b. 


.1 


bb) 


5 


# 


BR( 


~g 


-> 


~b. 


.1* 


b ) 


-5 


# 


BR( 


~g 


-> 


~b_ 


.2 


bb) 


5 


# 


BR( 


~g 


-> 


~b. 


2* 


b ) 


-6 


# 


BR( 


~g 


-> 


~t. 


.1 


tb) 


6 


# 


BR( 


~g 


-> 


~t. 


.1* 


t ) 


-6 


# 


BR( 


~g 


-> 


~t. 


.2 


tb) 


6 


# 


BR( 


~g 


-> 


~t. 


.2* 


t ) 



# stopl decays 
ID2 



6 


# BR(~t. 


.1 


-> 


~chi. 


.10 


t 


) 


6 


# BR(~t. 


.1 


-> 


~chi. 


.20 


t 


) 


6 


# BR(~t. 


.1 


-> 


~chi. 


.30 


t 


) 


6 


# BR(~t. 


.1 


-> 


~chi. 


.40 


t 


) 


5 


# BR(~t. 


.1 


-> 


~chi. 


.1+ 


b 


) 


5 


# BR(~t. 


.1 


-> 


~chi. 


.2+ 


b 


) 


6 


# BR(~t. 


.1 


-> 


~g 




t 


) 



24 






OOOOOOOOE+00 




2 


1000005 




37 


# BR(~t. 


.1 


-> 


~b_l 




h+: 





OOOOOOOOE+00 




2 


2000005 




37 


# BR(~t. 


_1 


-> 


~b_2 




H+, 





OOOOOOOOE+00 




2 


1000005 




24 


# BR(~t. 


.1 


-> 


~b_l 




W+ y 





OOOOOOOOE+00 




2 


2000005 




24 


# BR(~t. 


.1 


-> 


~b_2 




T T . > 

W+ y 


# 
# 


PDG 






Width 


















DECAY 


2000006 


7 


07039620E+00 


# 


stop2 decays 












# 


BR 




Tim A 

JNDA 


1JJ1 




lL)z 














2 


63043489E-02 




o 
z 


1 AAAAOO 

lUUUUzz 







4+ TDD ( ~+- 
if Dill, t. 


_2 


-> 


~chi. 


.10 


t ] 


9 


04476605E-02 




z 


-1 AAAAOQ 

lUUUUzo 




c 
D 


4f DD t ~ +- 


_2 


-> 


~chi. 


.20 


t ) 


4 


42222844E-02 




o 
z 


1 A A A AO C 

lvvvvzo 







4+ DD ( ~4- 
if Dlil, t. 


_2 


-> 


~chi. 


.30 


t ) 


1 


96823757E-01 




2 


1000035 




6 


# BR(~t. 


_2 


-> 


~chi. 


.40 


t ) 


2 


29818044E-01 




2 


1000024 




5 


# BR(~t. 


_2 


-> 


~chi. 


.1+ 


b ] 


2 


01693594E-01 




o 
z 


-\ AAAAQV 

lUUUUo 1 




b 


4+ DD ( ~ +- 

if dKI, t. 


_2 


-> 


~chi. 


.2+ 


b ) 





00000000E+00 




2 


1000021 




D 


4f DD t ~ +- 
if Dli\ L_ 


_2 


-> 


~g 




t ] 


3 


81187184E-02 




2 


1000006 




zb 


4+ DD ( ~+- 

ff brtA t. 


_2 


-> 


~t_l 




h ) 





00000000E+00 




2 


1000006 




ob 


44 DD ^ 

tt dKI, t. 


_2 


-> 


~t_l 




H ) 





00000000E+00 




2 


1000006 




OD 


44 DD ( ~ +- 
if b-Kv. t. 


_2 


-> 


~t_l 




A ] 





00000000E+00 




2 


1000005 




"27 


44 DD ( ~4- 


_2 


-> 


~b_l 




h+: 





00000000E+00 




2 


2000005 




Q7 


44 DD ( ~4- 


_2 


-> 


~b_2 




h+: 


1 


72571593E-01 




2 


1000006 




0*3 
ZO 


44 DD ^ ~4- 


_2 


-> 


~t_l 




Z ) 





00000000E+00 




2 


1000005 




O/l 

z4 


44 DD ( ~+~ 


_2 


-> 


~b_l 









00000000E+00 




2 


2000005 






44 DD ( ~4- 


_2 


-> 


~b_2 




W+^ 


# 
# 


PDG 






Width 


















DECAY 


1000005 


3 


79216587E+00 


# 


sbottoml 


decays 












# 


BR 




NDA 


ID1 




ID2 














4 


61784398E-02 




2 


1000022 




5 


# BR(~b. 


.1 


-> 


~chi. 


.10 


b ) 


3 


44388596E-01 




2 


1000023 




5 


# BR(~b. 


.1 


-> 


~chi. 


.20 


b ) 


5 


08897163E-03 




2 


1000025 




5 


# BR(~b. 


.1 


-> 


~chi. 


.30 


b ) 


1 


06407337E-02 




2 


1000035 




5 


# BR(~b. 


_1 


-> 


~chi. 


.40 


b ) 


4 


49199440E-01 




2 


-1000024 




6 


# BR(~b. 


.1 


-> 


~chi. 


.1- 


t ] 





00000000E+00 




2 


-1000037 




6 


# BR(~b. 


_1 


-> 


~chi. 


.2- 


t ] 





00000000E+00 




2 


1000021 




5 


# BR(~b. 


.1 


-> 


~g 




b ; 





00000000E+00 




2 


1000006 




-37 


# BR(~b. 


.1 


-> 


~t_l 




h-; 





00000000E+00 




2 


2000006 




-37 


# BR(~b. 


.1 


-> 


~t_2 




h-: 


1 


44503819E-01 




2 


1000006 




-24 


# BR(~b. 


.1 


-> 


~t_l 




w-: 





00000000E+00 




2 


2000006 




-24 


# BR(~b. 


.1 


-> 


~t_2 




w-: 



# 

# PDG Width 

DECAY 2000005 9 . 30237022E-01 # sbottom2 decays 



25 





BR 


NDA 


ID1 


ID2 
















2 


18042454E-01 


2 


1000022 


5 


# 


BR(~b_ 


_2 


-> 


~chi_ 


.10 


b ; 


1 


63625100E-01 


2 


1000023 


5 


# 


BR(~b_ 


_2 


-> 


~chi_ 


.20 


b ; 


4 


87739508E-02 


2 


1000025 


5 


# 


BR(~b_ 


_2 


-> 


~chi_ 


.30 


b ; 


7 


28430988E-02 


2 


1000035 


5 


# 


BR(~b. 


_2 


-> 


~chi_ 


.40 


b ; 


2 


17719337E-01 


2 


-1000024 


6 


# 


BR(~b. 


_2 


-> 


~chi. 


.1- 


t ; 





OOOOOOOOE+00 


2 


-1000037 


6 


# 


BR(~b. 


_2 


-> 


~chi. 


.2- 


t ; 





OOOOOOOOE+00 


2 


1000021 


5 


# 


BR(~b. 


_2 


-> 


~g 




b ; 





OOOOOOOOE+00 


2 


1000005 


25 


# 


BR(~b. 


_2 


-> 


~b_l 




h ; 





OOOOOOOOE+00 


2 


1000005 


35 


# 


BR(~b. 


_2 


-> 


~b_l 




h ; 





OOOOOOOOE+00 


2 


1000005 


36 


# 


BR(~b. 


_2 


-> 


~b_l 




a ; 





OOOOOOOOE+00 


2 


1000006 


-37 


# 


BR(~b. 


_2 


-> 


~t_l 




h-; 





OOOOOOOOE+00 


2 


2000006 


-37 


# 


BR(~b. 


_2 


-> 


~t_2 




h-; 





OOOOOOOOE+00 


2 


1000005 


23 


# 


BR(~b. 


.2 


-> 


~b_l 




z ; 


2 


78996060E-01 


2 


1000006 


-24 


# 


BR(~b. 


_2 


-> 


~t_l 




w-; 





OOOOOOOOE+00 


2 


2000006 


-24 


# 


BR(~b. 


_2 


-> 


~t_2 




w-: 



# 

# PDG Width 

DECAY 1000002 5 . 56209623E+00 # sup.L decays 





BR 


NDA 


ID1 


ID2 














5 


04724246E-03 


2 


1000022 


2 


# BR( 


"u_ 


_L 


-> 


~chi_ 


.10 u) 


3 


17395881E-01 


2 


1000023 


2 


# BR( 


"u_ 


_L 


-> 


~chi_ 


.20 u) 


9 


89628847E-04 


2 


1000025 


2 


# BR( 


"u_ 


.L 


-> 


~chi. 


.30 u) 


1 


15831154E-02 


2 


1000035 


2 


# BR( 


"u_ 


_L 


-> 


~chi. 


.40 u) 


6 


50599773E-01 


2 


1000024 


1 


# BR( 


"u_ 


_L 


-> 


~chi. 


.1+ d) 


1 


43843593E-02 


2 


1000037 


1 


# BR( 


"u_ 


_L 


-> 


~chi_ 


.2+ d) 





00000000E+00 


2 


1000021 


2 


# BR( 


"u_ 


_L 


-> 


~g 


u) 



# 

# PDG Width 

DECAY 2000002 1 . 05910452E+00 # sup.R decays 





BR 


NDA 


ID1 


ID2 














9 


83357095E-01 


2 


1000022 


2 


# brC 


"u_ 


_R 


-> 


~chi. 


.10 u) 


1 


10363000E-02 


2 


1000023 


2 


# BRC 


"u_ 


_R 


-> 


~chi. 


.20 u) 


1 


32287897E-03 


2 


1000025 


2 


# BRC 


"u_ 


_R 


-> 


~chi. 


.30 u) 


4 


28372575E-03 


2 


1000035 


2 


# BRC 


"u_ 


_R 


-> 


~chi. 


.40 u) 





00000000E+00 


2 


1000024 


1 


# BRC 


~u_ 


_R 


-> 


~chi. 


.1+ d) 





00000000E+00 


2 


1000037 


1 


# BRC 


"u_ 


_R 


-> 


~chi. 


.2+ d) 





00000000E+00 


2 


1000021 


2 


# BRC 


"u_ 


_R 


-> 


~g 


u) 



# 

# PDG Width 

DECAY 1000001 5 . 35470938E+00 # sdown.L decays 
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# BR 
2.35733116E-02 
3.07006812E-01 
1.72033892E-03 
1.64027859E-02 
6.04624253E-01 
4.66724989E-02 
. OOOOOOOOE+00 

# 

# PDG 
DECAY 2000001 

# BR 
9.83395634E-01 
1.10304127E-02 
1.31629502E-03 
4.25765842E-03 
. 00000000E+00 
. 00000000E+00 
. 00000000E+00 

# 

# PDG 
DECAY 1000004 

# BR 

5 . 04724246E-03 
3.17395881E-01 
9.89628847E-04 
1.15831154E-02 
6.50599773E-01 
1.43843593E-02 
. 00000000E+00 

# 

# PDG 
DECAY 2000004 

# BR 
9.83357095E-01 
1.10363000E-02 
1.32287897E-03 
4.28372575E-03 
. 00000000E+00 
. 00000000E+00 
0. 00000000E+00 



NDA 


ID1 


2 


1000022 


2 


1000023 


2 


1000025 


2 


1000035 


2 


-1000024 


2 


-1000037 


2 


1000021 



Width 



2.64248238E-01 



NDA 


ID1 


2 


1000022 


2 


1000023 


2 


1000025 


2 


1000035 


2 


-1000024 


2 


-1000037 


2 


1000021 



Width 
5.56209623E+00 



NDA 


ID1 


2 


1000022 


2 


1000023 


2 


1000025 


2 


1000035 


2 


1000024 


2 


1000037 


2 


1000021 



Width 



1.05910452E+00 



NDA 


ID1 


2 


1000022 


2 


1000023 


2 


1000025 


2 


1000035 


2 


1000024 


2 


1000037 


2 


1000021 



ID2 

1 # BR(~d. 

1 # BR(~d. 

1 # BR(~d. 

1 # BR(~d. 

2 # BR(~d. 
2 # BR(~d. 
1 # BR(~d. 



# sdown_R decays 
ID2 

1 # BR(~d. 

1 # BR(~d. 

1 # BR(~d. 

1 # BR(~d. 

2 # BR(~d. 
2 # BR(~d. 
1 # BR(~d. 



# scharm_L decays 
ID2 

4 # BR(~c 

4 # BR(~c 

4 # BR(~c 

4 # BR(~c 

3 # BR(~c 

3 # BR(~c 

4 # BR(~c 



# scharm_R decays 
ID2 

4 # BR(~c 

4 # BR(~c 

4 # BR(~c 

4 # BR(~c 

3 # BR(~c 

3 # BR(~c 

4 # BR(~c 



L 


-> 


~chi. 


.10 


d) 


L 


-> 


~chi. 


.20 


d) 


L 


-> 


~chi. 


.30 


d) 


L 


-> 


~chi. 


.40 


d) 


L 


-> 


~chi. 


.1- 


u) 


L 


-> 


~chi. 


.2- 


u) 


L 


-> 


~g 




d) 



R 


-> 


~chi. 


.10 


d) 


R 


-> 


~chi. 


.20 


d) 


R 


-> 


~chi. 


.30 


d) 


R 


-> 


~chi. 


.40 


d) 


R 


-> 


~chi. 


.1- 


u) 


R 


-> 


~chi. 


.2- 


u) 


R 


-> 


~g 




d) 



L 


-> 


~chi. 


.10 


c) 


L 


-> 


~chi. 


.20 


c) 


L 


-> 


~chi. 


.30 


c) 


L 


-> 


~chi. 


.40 


c) 


L 


-> 


~chi. 


.1+ 


s) 


L 


-> 


~chi. 


.2+ 


s) 


L 


-> 


~g 




c) 



R 


-> 


~chi. 


.10 


c) 


R 


-> 


~chi. 


.20 


c) 


R 


-> 


~chi. 


.30 


c) 


R 


-> 


~chi. 


.40 


c) 


R 


-> 


~chi. 


.1+ 


s) 


R 


-> 


~chi. 


.2+ 


s) 


R 


-> 


~g 




c) 
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# 
# 

DECAY 
# 



# 
# 

DECAY 
# 



# 
# 

DECAY 
# 



# 
# 

DECAY 
# 



PDG 
1000003 
BR 

2.35733116E-02 
3.07006812E-01 
1.72033892E-03 
1.64027859E-02 
6.04624253E-01 
4.66724989E-02 
. 00000000E+00 

PDG 
2000003 
BR 

9.83395634E-01 
1.10304127E-02 
1.31629502E-03 
4.25765842E-03 
. 00000000E+00 
. 00000000E+00 
. 00000000E+00 

PDG 
1000011 
BR 

4.19337361E-01 
2.11871941E-01 
. 00000000E+00 
. 00000000E+00 
3.68790698E-01 
0.00000000E+00 

PDG 
2000011 
BR 

1 . 00000000E+00 
. 00000000E+00 
. 00000000E+00 
. 00000000E+00 
0.00000000E+00 



Width 
5.35470938E+00 



# sstrange_L decays 



NDA 


ID1 


ID2 














2 


1000022 


3 


# BR(~s. 


_L 


-> 


~chi_ 


.10 


s) 


2 


1000023 


3 


# BR(~s. 


_L 


-> 


~chi_ 


.20 


s) 


2 


1000025 


3 


# BR(~s. 


_L 


-> 


~chi. 


.30 


s) 


2 


1000035 


3 


# BR(~s. 


.L 


-> 


~chi. 


.40 


s) 


2 


-1000024 


4 


# BR(~s. 


_L 


-> 


~chi. 


.1- 


c) 


2 


-1000037 


4 


# BR(~s. 


_L 


-> 


~chi_ 


.2- 


c) 


2 


1000021 


3 


# BR(~s. 


.L 


-> 


~g 




s) 



Width 
2.64248238E-01 



# sstrange_R decays 



NDA 


ID1 


ID2 
















2 


1000022 


3 


# 


BR(~s. 


_R 


-> 


~chi. 


.10 


s) 


2 


1000023 


3 


# 


BR(~s. 


_R 


-> 


~chi. 


.20 


s) 


2 


1000025 


3 


# 


BR(~s. 


_R 


-> 


~chi. 


.30 


s) 


2 


1000035 


3 


# 


BR(~s. 


_R 


-> 


~chi. 


.40 


s) 


2 


-1000024 


4 


# 


BR(~s. 


_R 


-> 


~chi. 


.1- 


c) 


2 


-1000037 


4 


# 


BR(~s. 


.R 


-> 


~chi. 


2- 


c) 


2 


1000021 


3 


# 


BR(~s. 


.R 


-> 


~g 




s) 



Width 
2.53683595E-01 



NDA 


ID1 


2 


1000022 


2 


1000023 


2 


1000025 


2 


1000035 


2 


-1000024 


2 


-1000037 



# selectron_L decays 
ID2 



11 


# 


BR(~e_L 


-> 


~chi. 


.10 


e-) 


11 


# 


BR(~e_L 


-> 


~chi. 


.20 


e-) 


11 


# 


BR(~e_L 


-> 


~chi. 


.30 


e-) 


11 


# 


BR(~e_L 


-> 


~chi. 


.40 


e-) 


12 


# 


BR(~e_L 


-> 


~chi. 


.1- 


nu_e) 


12 


# 


BR(~e_L 


-> 


~chi. 


.2- 


nu_e) 



Width 
,93168017E-01 



NDA 
2 
2 
2 
2 
2 



ID1 
1000022 
1000023 
1000025 
1000035 
-1000024 



# selectron_R decays 
ID2 



11 


# 


BR(~e_R 


-> 


~chi_10 


e-) 


11 


# 


BR(~e_R 


-> 


~chi_20 


e-) 


11 


# 


BR(~e_R 


-> 


~chi_30 


e-) 


11 


# 


BR(~e_R 


-> 


~chi_40 


e-) 


12 


# 


BR(~e_R 


-> 


~chi_l- 


nu_e) 



28 



O.OOOOOOOOE+00 



# 
# 

DECAY 
# 



# 
# 

DECAY 
# 



# 
# 

DECAY 
# 



# 
# 

DECAY 
# 



PDG 
1000013 
BR 

19337361E-01 
.11871941E-01 
OOOOOOOOE+00 
00000000E+00 
68790698E-01 
OOOOOOOOE+00 

PDG 
2000013 
BR 

, OOOOOOOOE+00 
OOOOOOOOE+00 
, OOOOOOOOE+00 
. OOOOOOOOE+00 
, OOOOOOOOE+00 
OOOOOOOOE+00 

PDG 
1000015 
BR 

OOOOOOOOE+00 
OOOOOOOOE+00 
OOOOOOOOE+00 

. OOOOOOOOE+00 
OOOOOOOOE+00 
OOOOOOOOE+00 

. OOOOOOOOE+00 
OOOOOOOOE+00 

PDG 
2000015 
BR 

63906167E-01 
96659153E-01 
. OOOOOOOOE+00 
OOOOOOOOE+00 



2 -1000037 

Width 
2.53683595E-01 



12 # BR(~e_R -> ~chi_2- nu_e) 



NDA 
2 
2 
2 
2 
2 
2 



ID1 
1000022 
1000023 
1000025 
1000035 
-1000024 
-1000037 



# smuon_L decays 
ID2 

13 
13 
13 
13 
14 
14 



# BR(~mu_L -> 

# BR(~mu_L -> 

# BR(~mu_L -> 

# BR(~mu_L -> 

# BR(~mu_L -> 

# BR(~mu_L -> 



'chi_10 mu-) 
'chi_20 mu-) 
'chi_30 mu-) 
'chi_40 mu-) 
'chi_l- nu_mu) 
'chi_2- nu_mu) 



Width 
1.93168017E-01 



# smuon_R decays 



NDA 


ID1 


ID2 
















2 


1000022 


13 


# BRC 


"mu_R 


-> 


~chi. 


.10 


mu- 


■) 


2 


1000023 


13 


# BRC 


"mu_R 


-> 


~chi. 


.20 


mu- 


-) 


2 


1000025 


13 


# BRC 


"mu_R 


-> 


~chi. 


.30 


mu- 


-) 


2 


1000035 


13 


# BRC 


"mu_R 


-> 


~chi. 


.40 


mu- 


-) 


2 


-1000024 


14 


# BRC 


"mu_R 


-> 


"chi. 


.1- 


nu. 


.mu) 


2 


-1000037 


14 


# BRC 


"mu_R 


-> 


"chi. 


.2- 


nu. 


.mu) 



Width 
1.34568207E-01 



# stau_l decays 



NDA 


ID1 


ID2 










2 


1000022 


15 


# BR(~tau. 


.1 


-> 


~chi_10 


2 


1000023 


15 


# BR(~tau. 


.1 


-> 


~chi_20 


2 


1000025 


15 


# BR(~tau. 


.1 


-> 


~chi_30 


2 


1000035 


15 


# BR(~tau. 


.1 


-> 


~chi_40 


2 


-1000024 


16 


# BR(~tau. 


.1 


-> 


~chi_l- 


2 


-1000037 


16 


# BR(~tau. 


.1 


-> 


~chi_2- 


2 


1000016 


-37 


# BR(~tau. 


.1 


-> 


~nu_tauL 


2 


1000016 


-24 


# BR(~tau. 


.1 


-> 


~nu_tauL 



tau-) 

tau-) 

tau-) 

tau-) 

nu_tau) 

nu_tau) 



Width 

3.03123027E-01 # stau_2 decays 

NDA ID1 ID2 

2 1000022 15 # BR(~tau_2 -> ~chi_10 

2 1000023 15 # BR(~tau_2 -> ~chi_20 

2 1000025 15 # BR(~tau_2 -> ~chi_30 

2 1000035 15 # BR(~tau_2 -> ~chi_40 



tau-) 
tau-) 
tau-) 
tau-) 



29 



3.39434680E-01 
. OOOOOOOOE+00 
. OOOOOOOOE+00 
. OOOOOOOOE+00 
. OOOOOOOOE+00 
. OOOOOOOOE+00 
. OOOOOOOOE+00 
. OOOOOOOOE+00 

# 

# PDG 
DECAY 1000012 

# BR 
8.14895139E-01 
4.97403776E-02 
. 00000000E+00 
. 00000000E+00 
1.35364484E-01 
. 00000000E+00 

# 

# PDG 
DECAY 1000014 

# BR 
8.14895139E-01 
4.97403776E-02 
. 00000000E+00 
. 00000000E+00 
1.35364484E-01 
. 00000000E+00 

# 

# PDG 
DECAY 1000016 

# BR 
8.36412025E-01 
4.45980206E-02 
0. 00000000E+00 
. 00000000E+00 
1.18989955E-01 
. 00000000E+00 
. 00000000E+00 
. 00000000E+00 
0. 00000000E+00 



2 


-1000024 


2 


-1000037 


2 


1000016 


2 


1000016 


2 


1000015 


2 


1000015 


2 


1000015 


2 


1000015 



Width 
1.73484437E-01 # 
NDA ID1 
2 1000022 
2 1000023 
2 1000025 
2 1000035 
2 1000024 
2 1000037 

Width 
1.73484437E-01 # 
NDA ID1 
2 1000022 
2 1000023 
2 1000025 
2 1000035 
2 1000024 
2 1000037 



Width 
1.67039909E-01 # 



NDA 


ID1 


2 


1000022 


2 


1000023 


2 


1000025 


2 


1000035 


2 


1000024 


2 


1000037 


2 


1000015 


2 


2000015 


2 


1000015 



16 


# BR( tau_2 


-> 


~chi_l- 


nu_ 


16 


# BR( tau_2 


-> 


~chi_2- 


nu_ 


-37 


# BR( tau_2 


-> 


~nu_tauL 


H-) 


-24 


# BR(~tau_2 


-> 


~nu_tauL 


W-) 


25 


# BR(~tau_2 


-> 


~tau_l 


h) 


35 


# BR(~tau_2 


-> 


~tau_l 


H) 


36 


# BR(~tau_2 


-> 


~tau_l 


A) 


23 


# BR(~tau_2 


-> 


~tau_l 


Z) 



tau) 
tau) 



snu_eL decays 
ID2 



12 


# BR(" 


~nu_ 


.eL 


-> 


~chi. 


.10 


nu_e) 


12 


# BR(" 


~nu_ 


.eL 


-> 


~chi. 


.20 


nu_e) 


12 


# BRC 


~nu_ 


.eL 


-> 


~chi. 


.30 


nu_e) 


12 


# BRC 


~nu_ 


.eL 


-> 


~chi. 


.40 


nu_e) 


11 


# BRC 


"nu_ 


.eL 


-> 


~chi. 


.1+ 


e-) 


11 


# BRC 


~nu_ 


.eL 


-> 


~chi. 


.2+ 


e-) 



snu_muL decays 
ID2 



14 


# BRC 


"nu. 


_muL 


-> 


~chi. 


.10 


nu. 


_mu) 


14 


# BRC 


"nu. 


_muL 


-> 


~chi. 


.20 


nu. 


mu) 


14 


# BRC 


"nu. 


.muL 


-> 


~chi. 


.30 


nu. 


mu) 


14 


# BRC 


"nu. 


.muL 


-> 


~chi. 


.40 


nu. 


.mu) 


13 


# BRC 


"nu. 


.muL 


-> 


~chi. 


.1+ 


mu- 


-) 


13 


# BRC 


"nu. 


_muL 


-> 


~chi. 


.2+ 


mu- 


-) 



snu_tauL decays 
ID2 

16 # BR(~nu_tauL -> ~chi_10 nu.tau) 

16 # BR(~nu_tauL -> ~chi_20 nu.tau) 

16 # BR(~nu_tauL -> ~chi_30 nu.tau) 

16 # BR(~nu_tauL -> ~chi_40 nu.tau) 

15 # BR(~nu_tauL -> ~chi_l+ tau-) 

15 # BR(~nu_tauL -> ~chi_2+ tau-) 

-37 # BR(~nu_tauL -> ~tau_l+ H-) 

-37 # BR(~nu_tauL -> ~tau_2+ H-) 

-24 # BR(~nu_tauL -> ~tau_l+ W-) 



30 



O.OOOOOOOOE+00 

# 

# PDG 
DECAY 1000024 

# BR 

. 00000000E+00 
. 00000000E+00 
. 00000000E+00 
. 00000000E+00 
. 00000000E+00 
. 00000000E+00 
0.00000000E+00 
. 00000000E+00 
. 00000000E+00 
. 00000000E+00 
. 00000000E+00 
. 00000000E+00 
. 00000000E+00 
. 00000000E+00 
. 00000000E+00 
. 00000000E+00 
. 00000000E+00 
. 00000000E+00 
. 00000000E+00 
1 . 00000000E+00 
. 00000000E+00 
. 00000000E+00 
. 00000000E+00 
. 00000000E+00 
. 00000000E+00 
. 00000000E+00 
0.00000000E+00 
. 00000000E+00 
. 00000000E+00 

# 

# PDG 
DECAY 1000037 

# BR 

. 00000000E+00 
. 00000000E+00 
0.00000000E+00 



2 2000015 
Width 



1.15215450E-02 # 



NDA 


ID1 


2 


1000002 


2 


2000002 


2 


-1000001 


2 


-2000001 


2 


1000004 


2 


2000004 


2 


-1000003 


2 


-2000003 


2 


1000006 


2 


2000006 


2 


-1000005 


2 


-2000005 


2 


1000012 


2 


1000014 


2 


1000016 


2 


-1000011 


2 


-2000011 


2 


-1000013 


2 


-2000013 


2 


-1000015 


2 


-2000015 


2 


1000022 


2 


1000023 


2 


1000025 


2 


1000035 


2 


1000022 


2 


1000023 


2 


1000025 


2 


1000035 



Width 
2.51070848E+00 # 
NDA ID1 
2 1000002 
2 2000002 
2 -1000001 



-24 # BR(~nu_tauL 



charginol+ decays 
ID2 



-1 


# 


BR( 


~chi_ 


1+ 


-1 


# 


BR( 


~chi_ 


1+ 


2 


# 


BR( 


~chi_ 


1+ 


2 


# 


BR( 


~chi_ 


1+ 


-3 


# 


BR( 


~chi_ 


1+ 


-3 


# 


BR( 


~chi_ 


1+ 


4 


# 


BR( 


~chi_ 


1+ 


4 


# 


BR( 


~chi_ 


1+ 


-5 


# 


BR( 


~chi_ 


1+ 


-5 


# 


BR( 


~chi_ 


1+ 


6 


# 


BR( 


~chi_ 


1+ 


6 


# 


BR( 


~chi_ 


1+ 


-11 


# 


BR( 


~chi_ 


1+ 


-13 


# 


BR( 


~chi_ 


1+ 


-15 


# 


BR( 


~chi_ 


1+ 


12 


# 


BR( 


~chi_ 


1+ 


12 


# 


BR( 


~chi_ 


1+ 


14 


# 


BR( 


~chi_ 


1+ 


14 


# 


BR( 


~chi_ 


1+ 


16 


# 


BR( 


~chi_ 


1+ 


16 


# 


BR( 


~chi_ 


1+ 


24 


# 


BR( 


~chi_ 


1+ 


24 


# 


BR( 


~chi_ 


1+ 


24 


# 


BR( 


~chi_ 


1+ 


24 


# 


BR( 


~chi_ 


1+ 


37 


# 


BR( 


~chi_ 


1+ 


37 


# 


BR( 


~chi_ 


1+ 


37 


# 


BR( 


~chi_ 


1+ 


37 


# 


BR( 


~chi_ 


1+ 



chargino2+ decays 
ID2 

-1 # BR(~chi_2+ 
-1 # BR(~chi_2+ 
2 # BR(~chi_2+ 



-> ~tau_2+ W-) 



-> ~u_L db) 

-> ~u_R db) 

-> ~d_L* u ) 

-> ~d_R* u ) 

-> ~c_L sb) 

-> ~c_R sb) 

-> ~s_L* c ) 

-> ~s_R* c ) 

-> ~t_l bb) 

-> ~t_2 bb) 

-> ~b_l* t ) 

-> ~b_2* t ) 

-> ~nu_eL e+ ) 

-> ~nu_muL mu+ ) 

-> ~nu_taul tau+) 

-> ~e_L+ nu_e) 

-> ~e_R+ nu_e) 

-> ~mu_L+ nu_mu) 

-> ~mu_R+ nu_mu) 

-> ~tau_l+ nu_tau) 

-> ~tau_2+ nu_tau) 

-> ~chi_10 W+) 

-> ~chi_20 W+) 

-> ~chi_30 W+) 

-> ~chi_40 W+) 

-> ~chi_10 H+) 

-> ~chi_20 H+) 

-> ~chi_30 H+) 

-> ~chi_40 H+) 



-> 


~u_L 


db) 


-> 


~u_R 


db) 


-> 


~d_L* 


u ) 






OOOOOOOOE+00 


2 


-2000001 


2 


# 


BR( 


~chi. 


.2+ 


-> 


~d_R* 


u ) 





OOOOOOOOE+00 


2 


1000004 


-3 


# 


BR( 


~chi_ 


.2+ 


-> 


~c_L 


sb) 





OOOOOOOOE+00 


2 


2000004 


-3 


# 


BR( 


~chi. 


.2+ 


-> 


~c_R 


sb) 





OOOOOOOOE+00 


2 


-1000003 


4 


# 


BR( 


~chi. 


.2+ 


-> 


~s_L* 


c ) 





OOOOOOOOE+00 


2 


-2000003 


4 


# 


BR( 


~chi. 


.2+ 


-> 


~s_R* 


c ) 





OOOOOOOOE+00 


2 


1000006 


-5 


# 


BR( 


~chi. 


.2+ 


-> 


~t_l 


bb) 





OOOOOOOOE+00 


2 


2000006 


-5 


# 


BR( 


~chi_ 


.2+ 


-> 


~t_2 


bb) 





OOOOOOOOE+00 


2 


-1000005 


6 


# 


BR( 


~chi. 


.2+ 


-> 


~b_l* 


t ) 





OOOOOOOOE+00 


2 


-2000005 


6 


# 


BR( 


~chi. 


.2+ 


-> 


~b_2* 


t ) 


1 


96925255E-02 


2 


1000012 


-11 


# 


BR( 


~chi. 


.2+ 


-> 


~nu_eL 


e+ ) 


1 


96925255E-02 


2 


1000014 


-13 


# 


BR( 


~chi. 


.2+ 


-> 


~nu_muL 


mu+ ) 


2 


66937138E-02 


2 


1000016 


-15 


# 


BR( 


~chi_ 


.2+ 


-> 


~nu_taul 


tau+) 


5 


23600711E-02 


2 


-1000011 


12 


# 


BR( 


~chi. 


.2+ 


-> 


~e_L+ 


nu_e) 





OOOOOOOOE+00 


2 


-2000011 


12 


# 


BR( 


~chi. 


.2+ 


-> 


~e_R+ 


nu_e) 


5 


23600711E-02 


2 


-1000013 


14 


# 


BR( 


~chi. 


.2+ 


-> 


~mu_L+ 


nu_mu ) 





OOOOOOOOE+00 


2 


-2000013 


14 


# 


BR( 


~chi. 


.2+ 


-> 


~mu_R+ 


nu_mu ) 


2 


27029940E-04 


2 


-1000015 


16 


# 


BR( 


~chi_ 


.2+ 


-> 


~tau_l+ 


nu_tau) 


5 


67355363E-02 


2 


-2000015 


16 


# 


BR( 


~chi. 


.2+ 


-> 


~tau_2+ 


nu_tau) 


2 


42517968E-01 


2 


1000024 


23 


# 


BR( 


~chi. 


.2+ 


-> 


~chi_l+ 


z ) 


6 


41934286E-02 


2 


1000022 


24 


# 


BR( 


~chi. 


.2+ 


-> 


~chi_10 


w+) 


2 


95919427E-01 


2 


1000023 


24 


# 


BR( 


~chi. 


.2+ 


-> 


~chi_20 


w+) 





OOOOOOOOE+00 


2 


1000025 


24 


# 


BR( 


~chi_ 


.2+ 


-> 


~chi_30 


w+) 





OOOOOOOOE+00 


2 


1000035 


24 


# 


BR( 


~chi. 


.2+ 


-> 


~chi_40 


w+) 


1 


69607703E-01 


2 


1000024 


25 


# 


BR( 


~chi. 


.2+ 


-> 


~chi_l+ 


h ) 





OOOOOOOOE+00 


2 


1000024 


35 


# 


BR( 


~chi. 


.2+ 


-> 


~chi_l+ 


H ) 





OOOOOOOOE+00 


2 


1000024 


36 


# 


BR( 


~chi. 


.2+ 


-> 


~chi_l+ 


A ) 





OOOOOOOOE+00 


2 


1000022 


37 


# 


BR( 


~chi_ 


.2+ 


-> 


~chi_10 


H+) 





OOOOOOOOE+00 


2 


1000023 


37 


# 


BR( 


~chi. 


.2+ 


-> 


~chi_20 


H+) 





OOOOOOOOE+00 


2 


1000025 


37 


# 


BR( 


~chi. 


.2+ 


-> 


~chi_30 


H+) 





OOOOOOOOE+00 


2 


1000035 


37 


# 


BR( 


~chi. 


.2+ 


-> 


~chi_40 


H+) 



# 

# PDG Width 

DECAY 1000022 . 00000000E+00 # neutralinol decays 
# 

# PDG Width 

DECAY 1000023 1 . 64735255E-02 # neutralino2 decays 





BR 


NDA 


ID1 


ID2 





















OOOOOOOOE+00 


2 


1000022 


23 


# 


BR(~ 


chi. 


.20 


-> 


~chi. 


.10 


Z ) 





OOOOOOOOE+00 


2 


1000024 


-24 


# 


BR(~ 


chi. 


.20 


-> 


~chi. 


.1+ 


W-) 





OOOOOOOOE+00 


2 


-1000024 


24 


# 


BR(~ 


chi. 


.20 


-> 


~chi. 


.1- 


W+) 





OOOOOOOOE+00 


2 


1000037 


-24 


# 


BR(~ 


chi. 


.20 


-> 


~chi. 


.2+ 


W-) 



32 






OOOOOOOOE+00 


2 


-1000037 


24 


# 


BR( 


[~chi. 


.20 


-> 


~chi_2- 


W+) 





OOOOOOOOE+00 


2 


1000022 


25 


# 


BR( 


>~chi_ 


.20 


-> 


~chi_10 


h ) 





OOOOOOOOE+00 


2 


1000022 


35 


# 


BR( 


i~chi_ 


.20 


-> 


~chi_10 


H ) 





OOOOOOOOE+00 


2 


1000022 


36 


# 


BR( 


i~chi_ 


.20 


-> 


~chi_10 


A ) 





OOOOOOOOE+00 


2 


1000024 


-37 


# 


BR( 


>~chi_ 


.20 


-> 


~chi_l+ 


H-) 





OOOOOOOOE+00 


2 


-1000024 


37 


# 


BR( 


.~chi. 


.20 


-> 


~chi_l- 


H+) 





OOOOOOOOE+00 


2 


1000037 


-37 


# 


BR( 


i~chi_ 


.20 


-> 


~chi_2+ 


H-) 





OOOOOOOOE+00 


2 


-1000037 


37 


# 


BR( 


.~chi. 


.20 


-> 


~chi_2- 


H+) 





OOOOOOOOE+00 


2 


1000002 


-2 


# 


BR( 


;~chi_ 


.20 


-> 


~u_L 


ub) 





OOOOOOOOE+00 


2 


-1000002 


2 


# 


BR( 


'~chi_ 


.20 


-> 


~u_L* 


u ) 





OOOOOOOOE+00 


2 


2000002 


-2 


# 


BR( 


.~chi. 


.20 


-> 


~u_R 


ub) 





OOOOOOOOE+00 


2 


-2000002 


2 


# 


BR( 


i~chi_ 


.20 


-> 


~u_R* 


u ) 





OOOOOOOOE+00 


2 


1000001 


-1 


# 


BR( 


.~chi. 


.20 


-> 


~d_L 


db) 





OOOOOOOOE+00 


2 


-1000001 


1 


# 


BR( 


i~chi_ 


.20 


-> 


~d_L* 


d ) 





OOOOOOOOE+00 


2 


2000001 


-1 


# 


BR( 


>~chi_ 


.20 


-> 


~d_R 


db) 





OOOOOOOOE+00 


2 


-2000001 


1 


# 


BR( 


.~chi. 


.20 


-> 


~d_R* 


d ) 





OOOOOOOOE+00 


2 


1000004 


-4 


# 


BR( 


>~chi_ 


.20 


-> 


~c_L 


cb) 





OOOOOOOOE+00 


2 


-1000004 


4 


# 


BR( 


>~chi_ 


.20 


-> 


~c_L* 


c ) 





OOOOOOOOE+00 


2 


2000004 


-4 


# 


BR( 


[~chi. 


.20 


-> 


~c_R 


cb) 





OOOOOOOOE+00 


2 


-2000004 


4 


# 


BR( 


'~chi_ 


.20 


-> 


~c_R* 


c ) 





OOOOOOOOE+00 


2 


1000003 


-3 


# 


BR( 


[~chi. 


.20 


-> 


~s_L 


sb) 





OOOOOOOOE+00 


2 


-1000003 


3 


# 


BR( 


i~chi_ 


.20 


-> 


~s_L* 


s ) 





OOOOOOOOE+00 


2 


2000003 


-3 


# 


BR( 


'~chi_ 


.20 


-> 


~s_R 


sb) 





OOOOOOOOE+00 


2 


-2000003 


3 


# 


BR( 


[~chi. 


.20 


-> 


~s_R* 


s ) 





OOOOOOOOE+00 


2 


1000006 


-6 


# 


BR( 


[~chi. 


.20 


-> 


~t_l 


tb) 





OOOOOOOOE+00 


2 


-1000006 


6 


# 


BR( 


[~chi. 


.20 


-> 


~t_l* 


t ) 





OOOOOOOOE+00 


2 


2000006 


-6 


# 


BR( 


i~chi_ 


.20 


-> 


~t_2 


tb) 





OOOOOOOOE+00 


2 


-2000006 


6 


# 


BR( 


[~chi. 


.20 


-> 


~t_2* 


t ) 





OOOOOOOOE+00 


2 


1000005 


-5 


# 


BR( 


;~chi_ 


.20 


-> 


~b_l 


bb) 





OOOOOOOOE+00 


2 


-1000005 


5 


# 


BR( 


i~chi_ 


.20 


-> 


~b_l* 


b ) 





OOOOOOOOE+00 


2 


2000005 


-5 


# 


BR( 


[~chi. 
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4.4 Some points about the program 

Our results for some representative points of the MSSM parameter space have been care- 
fully cross-checked against other existing codes in mSUGRA-type models [in the case of 
the AMSB and GMSB models, no very detailed comparisons have been made]. In our 
comparison with the programs Isajet and SPHEN0, when the sparticle and Higgs spec- 
trum, as well as the various input and soft SUSY-breaking parameters are forced to be the 
same in both codes, we obtain in general a rather good agreement for the partial widths 
of the main two-body decay modes. We have also verified that the loop and some three- 
body decays of the top squark agree qualitatively with those included in, respectively, the 
codes Isajet and SPHENO 6 . 

However, since there are, as is well known, some differences [of the order of a few 
per cent] between the outputs of the various RGE codes [see Ref. [44] for a discussion], 
these discrepancies can lead to a completely different phenomenology. Indeed, some decay 
channels can be either absent or present in some codes when the masses of the decaying 
and daughter particles are close to each other. For some multibody decays and for the 
QCD corrections, no comparison has been made since they are absent from these codes. 

The program is under rapid development and we plan to make several upgrades in a 
near future. A brief and non-exhaustive list of points that will be implemented in the 
next releases of the program includes: 

(i) the finalization of the inclusion of the decays of the top quark, i.e. the implementa- 
tion of the QCD corrections; 

(ii) the link with the routine Hmasses for a more precise and consistent [with the pro- 
gram] calculation of the Higgs boson masses; 
6 We thank Werner Porod for his gracious help in performing these detailed comparisons. 
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(Hi) the inclusion of additional higher order decay modes, such as the bottom squark 
three-body decay widths and some three-body decay modes, which can be important 
in GMSB models; 

(iv) the inclusion of the finite widths of the propagators of the exchanged particles in 
the multibody decay modes, to make smooth transitions between the multibody and 
two-body decay modes; 

(v) at some point, the implementation of some important electroweak radiative correc- 
tions [in particular for decays involving top and bottom squarks]. 

As mentioned earlier, a version where the program SDECAY is linked not only to the 
program SuSpect but also to the program HDECAY for the calculation of the decay widths 
and branching ratios of the MSSM Higgs bosons is in preparation. It would allow for a 
complete description of the properties of SUSY particles and MSSM Higgs bosons, prior 
to their production at colliders. 

Finally, a web page devoted to the SDECAY program can be found at the http address: 
http : //people . web . psi . ch/muehlleitner/SDECAY 
It contains all the information that is needed on the program: 

- downloading directly the various files of the program; 

- short explanations of the code and how to run it; 

- obtaining the complete "users manual" in post-script or PDF form; 

- a regularly updated list of important changes/corrected bugs in the code; 

- a mailing list to which one can subscribe to be automatically advised about future 
SDECAY updates or eventual corrections. 

5. Conclusions 

We have presented the Fortran code SDECAY, which evaluates the decay widths and branch- 
ing ratios of the supersymmetric particles in the MSSM. It includes not only all the pos- 
sible tree-level decays into two-body final states, but also the various three-body modes 
for charginos, neutralinos, gluinos and top squarks, the four-body decays of the lightest 
top squark and the loop-induced decays of the lightest top squark, the next-to-lightest 
neutralino and the gluino. In addition, the QCD corrections to the tree-level two-body 
decays, which involve strongly interacting particles and the dominant electroweak correc- 
tions due to the running of the gauge and fermion Yukawa couplings are incorporated. 
Furthermore, we have included the decays of the NLSP in GMSB models and the standard 
and SUSY decay modes of the heavy top quark. 

The program uses the SuSpect code for the calculation of the spectrum and the various 
soft SUSY-breaking parameters, but it can be easily linked to any other RGE code. It 
can be also linked to the program HDECAY for the decays of the MSSM Higgs bosons, to 
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provide a complete picture of the new particles predicted in the MSSM, except for their 
production properties. The latter are dealt with by the Monte Carlo event generators 
with which SDECAY can be easily linked since, in particular, it generates an output a la 
SUSY Les Houches Accord. The program is user-friendly, flexible for the choice of options 
and approximations, and quite fast. It therefore allows for a rather accurate, reliable and 
efficient study of the phenomenology of the MSSM superparticles, including the possibility 
of large scans of the parameter space. 

The program is under rapid development and will be maintained regularly to include 
upgrades, improvements and potentially, corrections. Any suggestion, comment or com- 
plaint from the potential users will be welcome. 
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